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Abstract 

Correction factor methodologies have been 
developed which use steady experimental or analytical 
pressure or force data to correct steady and unsteady 
aerodynamic calculations. Three methods of calculating 
correction factors have been developed to match steady 
surface pressure distributions, to match airfoil section 
forces and moments, and to match total forces and 
moments. Data for a rectangular supercritical wing that 
was previously tested in the NASA Langley Research 
Center Transonic Dynamics Tunnel have been used to 
determine correction factors to match surface pressure 
distributions for a range of Mach numbers. These 
correction factors have also been applied to unsteady 
aerodynamic calculations and comparisons have been made 
with oscillatory experimental data for a range of reduced 
frequencies at several Mach numbers. 


Nomenclature 
c reference chord 

ACp (X derivative of lining pressure coefficient with 

respect to angle of attack 
Cg section lift coefficient 

C L total lift coefficient 

C m section pitching moment coefficient 

C M total pitching moment coefficient 

E % maximum allowable error + 100 

F section property or total force 

j objective function involving correction factors 

J objective function involving forces to be matched 

k reduced frequency 

nb number of boxes 

nf number of forces to be matched 

11 s number of strips 

x/c fraction of chord measured from leading edge 
Q arbitrary objective function weighting factor 

V velocity 

« angle of attack 

q fraction of semispan from wing root 

X correction factor 

o) frequency (rad/scc) 

Subscripts 

D value of desired force to be matched 
i,k integer indices 

tot total 


Superscripts 

a analytical 

c experiment 

( )' derivative of ( ) with respect to ^ 

Introduction 

This paper describes the development of correction 
factor methodologies which use steady experimental or 
analytical pressure or force data to correct steady and 
unsteady aerodynamic calculations. The motivation for 
this research is that methods are still needed to improve 
routine analytical calculations of steady and unsteady 
pressures and forces using high quality aerodynamic data. 
Correction factors are multipliers which are applied to 
aerodynamic downwashes or pressures in aerodynamic 
calculations to achieve a specified desired objective. For 
example, figure 1 shows a typical aerodynamic box layout 
required to calculate unsteady aerodynamic forces using the 
Doublet Lattice Method (DLM) [1], where the correction 
factors would be applied to the individual boxes. 

Modifying analytical aerodynamic calculations 
using experimental aerodynamic data is not new. Gicsing, 
Kalman and Rodden [2] developed a correction factor 
technique to modify Doublet Lattice aerodynamics by 
matching total experimental force data. They formulated 
the problem such that the total forces were matched 
exactly while the change in the original analytical pressure 
distribution was minimized 

Pitt [3] recently developed a correction factor 
technique which modified strip aerodynamic influence 
coefficients (AIC's) using data generated by the transonic 
full potential fluid dynamics program FL028. Steady 
state section lift and moment were calculated by 
integrating the steady state pressures output from FL028 
over the chord. The derivatives with respect to angle of 
attack and control surface deflections were calculated using 
finite differences between separate cases of FL028 where 
the geometry input had been modified to reflect the change 
in shape. Pitt used the correction factors derived from the 
steady data to modify unsteady Doublet Lattice 
calculations and his corrected unsteady aerodynamics 
resulted in improved transonic flutter predictions. 

Three methods to calculate correction factors have 
been developed. The first approach described in this paper 
is to require a match between analytical and experimental 
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surface pressure distributions, the second approach requires 
airfoil section characteristics to be matched and the third 
matches total forces or integrated pressures. The first and 
second approaches require interpolation of experimental 
pressure data from the measurement stations to the 
analytical box locations. For the first approach, the 
correction factors arc calculated as the ratio of 
experimental to analytical pressure coefficients (ACp a ) at 
discrete points on the wing. The second approach is 
similar blit the correction factors arc ratios of section 
properties, such as Cq (I , C n i a , etc. Optionally, 

optimization techniques can be used to determine 
correction factors which minimize section property errors 
and/or minimize the change in the analytical pressure 
distribution. The third approach uses optimization 
techniques to determine correction factors so that total 
forces, moments, or control derivatives arc matched. This 
approach extends the work of Gicsing, Kalman, and 
Rodden [2] by reformulating the optimization problem and 
introducing additional objective functions and constraints. 

This paper describes the three methods of correction 
factor calculation but only presents results for the pressure 
distribution matching approach. The paper is organized in 
three sections. The first describes the methodologies, the 
second describes results and the final section presents 
some concluding remarks. 


Correction Factors Methodology 
Matching pressures 

One approach to correction factor calculation is to 
compare the steady experimental and analytical pressure 
coefficients at the aerodynamic box locations. Pressure 
correction factors are then calculated for each individual 
box as the ratio of the experimental to analytical lifting 
pressure coefficient derivative ACp w as 


pressure data and therefore have been used exclusively in 
this paper. One additional feature of cubic splines is that 
the first derivative (slope) of the interpolated curve is 
obtained as the evaluation of an analytical expression 
rather than by a numerical difference process, thereby 
reducing numerical errors in derivative calculations. 

The interpolation procedures will be described using 
the Rectangular Supercritical Wing (RSW) of references 8 
and 9 as an example since this wing was also used to 
generate the numerical results to be presented later. 
Figure 2 shows the analytical box layout for the RSW 
with the pressure measurement stations superimposed. 
The experimental data must be interpolated from the 
measurement stations to the quarter-chord pressure points 
at the midspan location of each of the analytical boxes. 
Figure 3 illustrates the interpolation process. First, the 
measured experimental data arc spline fit chordwise at the 
experimental spanwisc measurement stations (figure 3a) 
and interpolated to the chord stations used in the analytical 
model. The chordwise interpolated data are then spline fit 
spanwise and interpolated to the analytical span stations 
(figure 3b). 

To avoid extrapolation of pressures in regions 
beyond the measurement stations of the RSW, the 
following techniques were used. At the trailing edge, 
additional "data points" were added which explicitly forced 
the pressure differential there to be zero. To avoid 
extrapolation on the inboard section, the pressure 
distribution was assumed to be symmetric from tip to tip 
for the purpose of determining the spline fits. Finally, to 
avoid problems at the outboard edge of the wing, the 

experimental data was divided by the factor -\/(l - T| 2 ) 
prior to spline fitting and interpolation. The actual values 
of the interpolated experimental pressures at each 
analytical span location were then recovered by 
multiplying the values of the in terpolated modified 

pressures by the factor \/( 1 - T| 2 ) . 


This method therefore requires that experimental 
pressure data be interpolated to the analytical aerodynamic 
box locations. Surface splines [4] and one-dimensional 
numerical splines [5,6] were considered as alternatives for 
the pressure interpolations. Briefly, surface splines arc 
based on the small deflection equation of an infinite plate 
and have a discontinuous second derivative at the input 
data points. One-dimensional cubic splines arc based on 
the small deflection equation of an infinite beam and have 
a continuous second derivative and a discontinuous third 
derivative at the data input points. Surface splines arc 
commonly used to interpolate clastic deflections of wings 
[7[, but in the author’s experience one-dimensional cubic 
splines have proved more reliable in the interpolation of 


The above described interpolations are performed for 
steady data sets at different angles of attack and these data 
are then spline fit as a function of angle of attack. The 
analytical first derivative feature of one-dimensional cubic 
splines is then used to obtain the desired experimental 
lifting pressure coefficient derivative ACp a at the 
aerodynamic box locations. 

The correction factors are just the ratios of the 
interpolated experimental and analytical pressure data at 
each box location (equation 1). This is illustrated in 
figure 4 which shows representative interpolated 
experimental pressure distributions and the uncorrcctcd 
analytical pressure distributions for the RSW. 

Correction factors can also be calculated which are 
applied to box downwashes instead of pressures if 
experimental downwash data are available. The 
experimental downwash for the RSW was not measured 
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directly but could be calculated if both the airfoil shape 
and the deformation of the wing in the tunnel were 
known. The wing deformation of the RSW in the tunnel 
was not measured. Therefore, for this paper, a pseudo- 
cxperimcnlal downwash distribution was determined using 
the analytical aerodynamic influence coefficient (AIC) 
matrix and the measured experimental pressure 
distribution. Equation 2 stales the relationship between 
the box downwashes w and the pressure distribution 
ACp (t in terms of the AIC matrix 

{ w } - I AIC) I" AC Pa J- (2) 

The downwash correction factors arc the ratios of 
the pseudo-experimental downwashes to the analytical 
downwashes as 


The same correction factor is applied to all the 
aerodynamic boxes along a chordwise strip at the span 
location for which the section characteristics are valid. 

In most eases however, it is best to match both 
Cj a and C m(X simultaneously. One way this can be 

accomplished is by using correction factors which vary 
linearly along the chord section as 

X(5) = X'<5-5) + X (6) 

where X' is the slope of the line and X is the value of the 
correction factor at the moment reference center. X’ and 
X arc determined by solving the system of linear 
equations 



(3) 


Pressure corrections have the effect of modifying 
only the pressure on the box to which they arc applied, 
whereas a downwash correction factor on one box affects 
the pressures on all other boxes to varying degrees 
through the AIC matrix. The effect of pressure correction 
factors on unsteady pressures is to modify the magnitude 
of the unsteady pressures. The only modification to phase 
is a shift of 180 degrees when the correction factor is 
negative. Downwash correction factors affect both the 
magnitude and phase of unsteady analytical pressure 
distributions. Applying either the pressure correction 
factors calculated from equation 1 or the downwash 
correction factors calculated by equation 3 achieve a 
complete match of steady analytical and experimental data. 

Section Charact eri stics 


A second approach to correction factor calculation 
involves matching one or more airfoil section properties, 
for example, lift and moment, (Cj a and C ma , 

respectively). If only one section property (ie. Cj a or 
Cm a ) needs to be matched, the pressure correction factor 

is simply the ratio of the experimental to the analytical 
section property for each airfoil section, for example 



(4) 



(7) 


where A is defined by 


[A] = 


<4(5) 


a 


a 


<4(5) J *Cp («($-*)"<* 

“ 5 = 0 a 


5=1 


( 8 ) 


Experimental section Cj a can be calculated either 

by integrating ACp a along the chord using Gauss- 

Legendre Quadrature [10] for example, or as the derivative 
of the section lift with respect to angle of attack. 


Alternatively, optimization techniques can be used 
to calculate correction factors which arc different for each 
of the boxes along a chordwise strip. The correction 
factors for each strip arc obtained as the design variables in 
an optimization problem where the objective is to 
minimize the square error in the desired section property 
(properties). For the fih span station, the objective 
function is 


J.= 

i 


nf 

I 





(9) 


or 


C. 


a 


Because this approach matches section properties 
and not pressure distributions, the correction factors 
calculated by this method may result in an unrealistic 
analytical pressure distribution. This effect can be 
alleviated by 1) adding an additional term to the objective 
(5) function for the idi strip which minimizes the difference 
in the pressure distribution as 


3 



(10) 


nh . 

V= I a-x k ) 2 

;iml/or 2) hy including upper and lower bounds on the 
correction factors as inequality constraints in the 
optimization. For ease 1), the total objective function to 
be minimized for each strip is 


where the summation is over the total number of forces to 
be matched (CL a , CM a .CL5 etc). As in the previous 

method, changes in the pressure distribution can be 
minimized by appending a term 

nb 2 

J- « 4 > 


= J + Q j 

lot. i v i J i 

i 


where Qj is an arbitrary weighting factor. 


to equation 13, where nb here is the total number of 
(11) boxes, or by implementing inequality constraints which 
provide upper and lower bounds on the correction factors. 
Additionally, inequality constraints on individual forces 
can be formulated as 


Instead of solving separate optimization problems 
for each strip, the correction factors for all the strips could 
be calculated simultaneously by summing the objective 
functions for each strip as 

IIS 

J, <>i = J| J t«t (12) 



(15) 


lulalfmcs 

A third approach to correction factor calculation is 
to match total forces, total moments, or integrated 
pressures (eg. Cl (X , Cm (X . etc.). Figure 5 shows a bar 

chart comparing analytical force coefficients calculated 
using the Doublet Lattice Method (DLM) and the box 
layout of Figure I with experimental force data from 
reference 1 1 for the F/A-18. The goal of the methodology 
is to find a set of correction factors which improves the 
correlation of the analytical force results with known force 
data. 

This methodology expands on the work of Gicsing, 
Kalman, and Koddcn|2j which uses a closed form solution 
of an optimization problem where total forces and 
moments arc matched exactly using equality constraints 
and the objective function minimizes the variation of the 
modified pressure distribution from the original analysis. 
Correction factors arc applied to either the aerodynamic 
box downwashes or pressures. The methodology described 
in this paper provides for more options in the formulation 
of the optimization problem, and uses numerical 
optimization rather than a closed form solution. This 
total force methodology is basically an extension of the 
section properties methodology described earlier. 

The primary objective function for the total force 
matching problem is 

j %!,('•■* - p '\) 2 <>» 


in order to improve the matching of specific forces. 

Numerical optimization to calculate the correction 
factors to match total forces or section properties can be 
achieved using numerical optimization techniques, such as 
those which are available in the Automated Design 
Synthesis (ADS)[12] nonlinear programming software 
package. 


Results 

The pressure distribution matching method for 
correction factor calculation has been implemented and 
results for this methodology will be presented in this 
paper. The results arc for a Rectangular Supercritical 
Wing (RSW) that was tested in the NASA Langley 
Research Center Transonic Dynamics Tunnel. Reference 
8 presents the geometric and structural properties of the 
RSW and reference 9 gives tabular listings of steady 
pressure data for a range of Mach numbers and angles of 
attack, and unsteady oscillatory experimental pressure data 
for a range of Mach numbers, reduced frequencies and 
mean angles of attack. Figure 2, mentioned earlier, shows 
the aerodynamic box layout used in the analysis to 
calculate the analytical pressures. Superimposed arc the 
locations at which experimental steady and unsteady 
pressure data were measured. Steady experimental data 
were interpolated from the measurement stations to the 
analytical locations as described in the methodology 
section of the paper in order for pressure and downwash 
correction factors to be calculated. Unsteady measured data 
were also interpolated to analytical locations in order that 
direct comparisons could be made between the 
experimental pressures and the uncorrcctcd and corrected 
analytical pressures. 
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Figure 6 shows a typical chordwisc distribution of 
steady upper surface, lower surface, and lifting pressures at 
Mach 0.4, showing the distribution of pressures on this 
supercritical wing. Note the high lift on the rear half of 
the airfoil where supercritical wings typically have high 
cumber. Figure 7 shows the effects of anglc-of-attack 
variation on the lifting pressure distribution, with 
chordwisc distributions at midspan shown in figure 7a and 
spanwisc distributions at midchord shown in figure 7b. 
Tlie derivatives of these lifting pressures with respect to 
angle of attack arc the quantities which arc to be matched 
in this analysis. Figure 4 shows a typical comparison of 
experimental and analytical steady ACp a distributions 

from which pressure correction factors are calculated. 
Note that the pressures arc undcrprcdictcd by the DLM 
toward the leading edge and reasonably well predicted 
toward the trailing edge. 

Figure 8 shows pressure correction factors Q(p) that 
were calculated from steady pressure data at a Mach 
number of 0.4 for the RSW. The chordwisc distributions 
at several span stations and the spanwisc distributions at 
several chord locations arc provided to give an overall 
view of the behavior of the correction factors over the 
surface. Note that the pressure correction factors towards 
the trailing edge of the supercritical airfoil increase. This 
is basically due to the low magnitudes of the pressures in 
this region, where the ratio of the experimental to 
analytical pressures can be large even though the actual 
pressure prediction error is small. The figures also show 
that the correction factors can be negative as occurcd 
toward the outboard edge of the wing at about the mid- to 
three-quarter chord location. This may be due to 
variations in the flow around the wing tip caused by the 
unconventional chordwisc pressure distribution over a 
supercritical wing, as illustrated in Figure 6. This 
unconventional flow around the wing tip is not 
represented by the Doublet Lattice Method. Figure 9 

shows the downwash correction factors (% w ) which were 

calculated to match the steady experimental data for the 
same test condition as figure 8. Applying either the 
pressure or downwash correction factors results in a match 
between the analytical and experimental steady pressure 
distributions. 

Unsteady analytical pressures were calculated using 
the Doublet Lattice Method to obtain oscillatory 
aerodynamics at reduced frequencies corresponding to 
points where unsteady wind tunnel pressure data were 
available. These calculations were conducted so that direct 
comparisons could be made between the unsteady 
experimental pressure data and corrected and uncorrcctcd 
analytical calculations. Figure 10 shows a comparison of 
unsteady experimental pressure data (magnitude and phase) 
and uncorrcctcd and corrected analytical data for two 
reduced frequencies at Mach 0.4. The purpose of this 
comparison was to determine whether the correction 
factors based on steady data will improve unsteady 
pressure calculations over a reduced frequency range. 


Figures 10a and 10b compare the magnitudes of the 
unsteady pressures at a span station of T|=.525 for reduced 
frequencies of 0.309 (10 Hertz oscillation), and 0.618 (20 
Hertz oscillation). The corresponding phases are shown in 
Figures 10c and lOd. Both downwash and pressure 
correction factors resulted in improved prediction of the 
pressure magnitudes on the leading edge of the airfoil, 
with mixed results farther aft where the magnitudes arc 
much smaller. The effect of downwash correction factors 
on phase is generally beneficial, with pwtrer results in the 
same regions as for the pressure magnitudes. 

Correction factors based on steady pressure data 
were calculated for four Mach numbers ranging from 
0.266 to 0.8. The pressure correction factors arc shown in 
Figure 1 1 and the downwash correction factors are shown 
in Figure 12. This analysis was conducted to determine if 
correction factors calculated to match pressure data at one 
Mach number might be applicable over a Mach number 
range. As shown, both sets of correction factors vary 
with Mach number, indicating that for the RSW, a single 
set of correction factors calculated at one Mach number is 
not adequate over a wide range of Mach numbers. 

Figure 13 shows the unsteady pressure magnitudes 
for Mach numbers 0.266, 0.7, and 0.8 for an oscillation 
frequency of 10 Hertz, and figure 10a shows the 
corresponding data for Mach 0.4. Figures 14 and 10c 
show the phase of the unsteady pressures for the same test 
conditions. It is noted that the test medium for Mach 
0.266 was air while the test medium for the other Mach 
numbers was freon. Although the reduced frequencies arc 
not constant as they should be to isolate a Mach number 
effect, the reduced frequencies for Mach numbers 0.266, 
0.7 and 0.8 do not vary widely, so that the trends in 
Figures 13 and 14 can be attributed mostly to Mach 
number effects. Both types of correction factors, 
downwash and pressure, result in improved leading edge 
matching of the pressure magnitudes, including the 
transonic(shock) effects at Mach 0.8. This result is 
similar to that shown by Pitl[3]. The downwash 
correction factors also provided some improvement in the 
phase distributions. 


Concluding Remarks 

Three methods of correction factor calculation have 
been developed to match surface pressure distributions, to 
match section properties, and to match total forces, 
moments and control derivatives. Pressure and downwash 
correction factors to match surface pressure distributions 
have been calculated based on steady experimental data for 
a range of Mach numbers for a rectangular supercritical 
wing. Uncorrected and corrected analytical pressure data 
have been calculated and compared with steady and 
unsteady experimental pressure data at a range of Mach 
numbers and reduced frequencies. Correction factors are 
shown to improve the analytical calculation of unsteady 
pressures and could be used to modify analytical pressure 
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calculations lo account for some nonlinear aerodynamic 
c Heels. 
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Figure 1 . Typical aerodynamic box layout for 
Doublet Lattice Method. 
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^ Wind-Tunnel Wall 

Wing Pitch Axis 
(.46c) 


Figure 2. Analytical box layout with experimental 
pressure measurement stations of 
Rectangular Supercritical Wing. 
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• Chordwise interpolated experimental data 
o S panwise interpolated data 



a)Chordwise Interpolation 


b)Spanwise Interpolation 


Figure 3. Interpolation of experimental lifting pressures to analytical locations. 


a) Chordwise pressure distribution it b) Spanwise Pressure Distribution at 

n =*.025 (near wing root) x/c= .025 (near leading edge) 



Figure 4. Comparison of steady experimental and analytical steady pressure 
distributions to calculate correction factors. 
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■ Experimental 
□ Analytical 



-.1 0 . -.1 


Figure 5. Comparison of F-18 total force codficients 
(M=0.8). 


M=.4, a =.2! degrees 



x/c 

f igure 6. Typical chordwise distribution of upper, 

lower and differential pressure distributions 
for RSW. 


ri=.525 





Figure 7. Effect of angle of attack on steady lifting 
pressure distribution of RSW (M=0.4). 




Figure 8. Pressure correction factors for RSW at M=0.4. 
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Figure 9. Downwash correction factors for RSW at M=0.4. 




Figure 10. Effect of oscillation frequency on unsteady ACp a . 
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1'jj.urc 13. Comparisons of uncorrectcd and corrected I4 - Comparisons of uncorrected and corrected 

magnitude of ACp a for Mach numbers phase ot ACp a tor Mach numbers 0.266- 
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